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(54) Title: METHOD FOR CALCULATING THE VARIANCE 
(57) Abstract 

This applicaton describes a method and device for measuring 
a variation value, where a value representative or a variation 
between a quantity Q and a reference quantity Qref is determined 
by calculating error samples xi and then determining the value R 
representative of the variation from (6) where each value of wj 
is recursively determined from the previous value wj-u Through 
this recursive determination, a nested loop is avoided, whereby the 
amount of necessary calculation is greatly reduced. 
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METHOD FOR CALCULATING THE VARIANCE 
[Field of the invention] 

The field of the invention relates to a method and device 
for measuring a value that characterizes a variation 
between a quantity and a reference quantity, where said 
variation is measured with respect to a parameter on which 
both the quantity and reference quantity depend. An example 
of this is the determination of a variation between a 
signal and a reference signal, both of which depend on 
time . 

[Background of the invention] 

A well known measure of variation is the variance which is 
defined as the mean value of the square of the difference 
between a sample vaiue and the mean value of the samples. 
The standard deviation is the square root of the variance. 

Specifically in connection with the characterization of 
timing signals, a measure of the time variation as a 
function of integration time is known, which is referred to 
as TDEV (time deviation) and is e.g. defined in ETSI 
European Telecommunication Standard ETS 3004 62-1, April 
1997 or ITU-T Recommendation G.810, Series G: Transmission 
Systems and Media, August 1996. A time deviation TDEV is 
defined as 



TDEV(nr 0 )= ^/^( x i*2n- 2 ^ + « + ^) ) < 1 > 
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where the angle brackets denote an ensemble average, and 



x 



10 



15 



with i = 1/ 2, . . . . N, are sample values of a time error 
function x(t) and these samples are taken at equidistant 
time intervals. The time dependent error function x(t) is 
defined as the difference between a clock generating time 
T(t) and a reference clock generating time Tref(t). Here, 
is therefore to be understood as the absolute (abstract) 
time, whereas T is a signal generated by a clock that 
depends on t and itself also represents a time. The N 
sample values are sampled at equal intervals xq, such that 

Xi=x(i * x 0 ) f i = 1/ 2 , N. x 0 is the sampling period 

and x = n • xq is the observation interval. 



20 The ensemble average relates to the observation interval in 
the sense that an average is taken over possible triplets 

x i+2n/ x i+n' x i • 



25 



The above mentioned documents contain an estimator formula 
for TDEV, which is : 



TDEV(nT 0 ) = 



i 



W-3/H-1 



6r?(N- 3«+l) % 



-i2 



(2) 



30 where n = 1, 2, ... the integer part of N/3. 



[Problem underlying the invention] 
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5 The problem with the above estimator formula is that a 
calculation results in two nested FOR-loops for the 
summation under the square root, and in three nested FOR- 
loops when calculating the values TDEV for all values of n, 
which is often required. With large values of N, which is 

10 typically the case, this results in a large calculation 

burden which either leads to long calculation times or to 
an increased amount of hardware for coping with the 
calculation burden in a reasonable amount of time. As an 
example, in communication systems using the above described 

15 TDEV as a control parameter for synchronization control, 

this leads to more circuitry, which makes the systems more 
complicated and more expensive. 



20 



25 



The above mentioned problem is not restricted to systems 
calculating the deviation for a clock signal with respect 
to a reference clock signal over time, but will occur in 
any system based on the above principle of calculating a 
value representative of a variation based on the double sum 



f 2 0) 

E( X i + 2n" 2X i-H, +X i) 
i=f,0) 



(3) , 



J 

where the lower and upper boundaries of the inner sum 
running over i are respective functions (f a and f 2 ) of the 
outer variable j . 

30 [Object of the invention] 

The present invention has the object of providing a better 
method of measuring a value characteristic of the variation 
between a quantity and a reference quantity, said value 
35 being based on the above described double sum. 
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5 [Summary of the invention] 

This object is solved by the methods and devices described 
in the independent claims appended to the present 
application. Advantageous embodiments are described in the 
10 dependent claims. 

The present invention greatly simplifies the determination 
of a value representative of the variation between a 
quantity and a reference quantity by employing a recursive 
15 formula, such that the two nested FOR- loops mentioned above 
can be avoided. 

More specifically, the present invention defines a first 
value wi as a function of the difference values Xj, which 
20 represent differences between the quantity and the 

reference quantity, and determines the value R, which 
corresponds to the above mentioned double sum, through 

R = w 2 + £ w ~ {6 ) 

7 = 2 

25 

where each consecutive value of Wj is not calculated by 
running through a total sum, but is recursively determined 
from the respectively previous value Wj_]_. n is a value 
between 1 and the integer part of N/3. 

30 

In this way the present invention achieves a simpler and 
less time consuming calculation of the value R, such that 
the measurement method of the present invention enables a 
simpler and cheaper hardware, without any decrease in 
35 efficiency. More specifically, the calculation of the 

double sum of the above mentioned equation (3) requires 
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5 floating point operations in a number depending on N , 

whereas the method of the present invention only requires a 
number of floating point operations depending on N. 

Therefore, the present invention decreases the required 
10 processing capacity and processing time, to thereby lead to 
better and/or cheaper hardware. 

According to a preferred embodiment of the present 
invention, the calculation of w x is done in accordance with 



15 



20 



Z(x k+ 2n- 2x k + n +X k)< 
k = \ 

and the recursive calculation of Wj in accordance with 

W j = Vi + x J-i+*> ~ 3 */-i+2„ + 3 *;-i+« ~ x j-\ ' 

According to another preferred embodiment of the present 
invention, the calculation of wx is performed by 
introducing a help variable y, where 



25 Yk = x k+2n " 2x n+k + x k 



(8) 



and 



W 1=XX (9) 
k=l 

30 

such that the recursive formula for wj becomes 



wj = Wj_i + Yj-l+n ~ 



(10) 
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5 

According to further preferred embodiments, the values R, 
which depend on n, can be used to calculate a variance and 
a deviation value, either for individual values of n, or 
for all values of n from 1 to the integer part of N/3. 

10 

As another preferred embodiment of the present invention, 
the measurement method is applied to a communication system 
requiring synchronization, in which the deviation value 
that is calculated relates to the variation between a clock 

15 signal and a reference clock signal. According to another 

preferred embodiment, the present invention is applied to a 
system in which a surface profile level is compared to a 
reference profile level, and samples are measured at equal 
location intervals, such that the measured characteristic 

20 variation value relates to the spatial variation of the 
surface level. 

As already mentioned above, the present invention can be 
applied to a system calculating a value characteristic of a 
25 variation on the basis of the above mentioned double sum 
shown in equation (3) . 

[Brief description of the figures] 

30 Fig. la and lb show a flow-chart illustrating a first 

embodiment of the present invention; 

show a flow-chart illustrating a second 
embodiment of the present invention; 
and 

Fig. 3a and 3b show a flow-chart illustrating a fourth 

embodiment of the present invention. 



Fig. 2a and 2b 

35 
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5 Fig> 4 is a diagram illustrating the sampling 

number N, the measurement period Ntq, 

the sampling interval nxo, and the 

basic interval xq- 

10 Fig. 5a and 5b show a flow-chart illustrating an 

embodiment that is a modification of 
the embodiment shown in Fig. 1. 

Fig. 6a and 6b show a flow-chart illustrating a third 

embodiment of the present invention. 

Fig. 7a and 7b show a flow-chart illustrating an 

embodiment of the present invention 
that is a modification of the 
20 embodiment shown in Fig. 6. 

[Detailed disclosure] 

Fig. 1 shows a flow-chart of a first embodiment of the 
25 present invention, where Fig. la shows the first part and 

Fig. lb the second part. As can be seen, in a first step SI 
a loop variable i is set to a starting value of 1. Then, in 
step S2, a quantity Q, which depends on a given parameter t 
is measured at equidistant intervals x 0 , together with a 

30 reference quantity Qref associated with Q, where Qref also 
depends on the parameter t. In this way, samples Q ± = QU 
t 0 ) and Qrefi = Qref(i • x 0 ) are measured. Also, the 
difference between Qi and Qrefi is determined to thereby 
supply error samples x^Qi-Qrefi- It is sel f -understood 

35 that the error function can also be defined as Qrefi " Qi • 
It should be noted that t represents any possible 
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5 parameter, such as time, location, etc. The interval x 0 is 
naturally measured in the same units as the parameter t. 

Then, N samples are determined by running through the loop 
constituted by steps S2, S3 and S4, where step S3 

10 determines if i = N, and if not, then S4 increments the 
loop variable i to i + 1. N is a parameter that 
characterizes the measurement period, as will be explained 
in more detail in connection with Fig. 4. With regard to 
the present invention, it is of no importance how N is set 

15 or acquired, i.e. with respect to the invention, N can be 
seen as an arbitrary, given parameter. 

It should be noted that the quantity Q can be any 
measurable quantity, just as the parameter t can be any 

20 parameter on which a quantity can depend. As an example, Q 
can be a clock signal T and the parameter t can be the 
actual time, so that x would be the difference between a 
clock signal T(t) and a reference clock signal Tref(t) at 
equal spacing sample points defined by the time interval 

25 x 0 . On the other hand, the quantity Q can equally well be a 
surface level value and the parameter referred to as t a 
one dimensional location along said surface, such that x 
represents error values between surface level values and 
reference surface level values at points of equal distance 

30 t 0 along a direction in said surface. 

After the samples Xi are determined, a second loop using a 
loop variable k determines a value W x . More specifically, 
in step S5 the loop variable k is set to 1 and the value W x 
35 is set to 0. Then, a loop consisting of steps S6, Si and S8 
determines a value in accordance with the above 
mentioned equation (4) which is 
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n 



10 



20 



25 



30 



(4) 

k=l 



This is done by replacing the value W x by W 1 +Xj c+ 2 n -2x k+n +x k 
in step S6, then determining if k = n in step S7 and 
incrementing k to k+1 in step S8 if the result in S7 is 
negative . 



It should also be noted that the value for n is acquired in 
step S5. As will be explained in more detail in connection 
15 with Fig. 4, n is a parameter that characterizes the 

observation interval. With regard to the present invention, 
n can be acquired or set in any desirable or suitable way. 
For example, n can be a stored value that is simply read by 
the system in step S5. Or the user may be prompted to enter 
a desired value. Or the system may have a sub-routine for 
determining a value of n. As can be seen, it is of no 
importance to the invention how and where n is set, as long 
as the value of n lies between 1 and the integer part of 
N/3. 



Then, in step S81 it is determined if n is equal to the 
integer part of N/3, and if this is the case, then the 
value R, which represents the above mentioned double sum 



J 



(3) 



is determined as R = . This is shown in step S9 (see 
Fig. lb) . 
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15 



20 



On the other hand, if in step S81 it is determined that n 
is smaller than the integer part of N/3, then the 
processing proceeds to step S10, in which again the above 
mentioned value of R is set to w x 2 but then processing 
continues by setting the loop variable j to 2 in step Sll 
and then entering the loop as S12, S13, S14 and S15 to 
thereby determine the value R as 



N-3n+\ 2 
R = M> 1 2 + ± Wj (6) 

j = 2 



where each value Wj is determined recursively from the 
previous value of w, i.e. Wj_i by 



= wj_x + 2Cj_i + 3 n - 3xj_i+2n + 3x j-l+n ~ (7) 



This is done by calculating the new value of wj in 
accordance with the above mentioned equation (7) in steps 
S12, then replacing the value of R by R + w j 2 in step S13, 
asking if the value j has reached N - 3n+l in step S14 and 
25 incrementing j to j+1 in step S15 if the result of step S14 
is negative. 

The processing after step S14 or S9 runs together again, 
i.e. the value of R is then employed in any suitable way, 
30 e.g. output to a control means that responds to the measure 
of variation between Q and Qref expressed by R, or the 
value R is processed further, e.g. a variance value VAR is 
determined as 



BNSOCCIO <WO 9963*51 A 1 I > 



/ 



10 



20 



WO 99/63451 PCI7EP99/038S4 

11 

_ R 
VAR " 6n 2 (Af->z + l) 



and/or a deviation value is determined as DEV = (VAR) 1 / 2 or 



DEV= " R 



(12) 



The above mentioned embodiment calculates one value 
representative of a variation between a quantity Q and a 
reference quantity Qref associated with Q in dependence on 
a parameter n, where n determines the observation interval 
15 in the form of n • x 0 • n can be assigned or have any value 
between 1 and the integer part of N/3, where N indicates 
the measurement period N • xq • 



As can be seen , the above embodiment does not contain any 
nested loops, i.e. loops that are part of a larger outer 
loop. In this way, the above described double sum 



2] BX H *-2XH 11 +X i ) 

/=/,(/) 



(3) 



25 is calculated only from single, independent loops, so that 
the amount of necessary calculation only depends linearly 
on the sample size, i.e. only depends on N, in contrast to 
a conventional solution of the double sum, which would lead 
to a number of calculations depending on N 2 . Therefore, the 

30 above described embodiment achieves a better measurement 
method, as the amount of calculation is considerably 
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5 reduced/ namely by a factor of N, which in turn allows a 
far simpler and/or more efficient hardware implementation. 

The present invention is by no means restricted to the 
specific embodiment described in Fig. la and lb, as a 

10 person skilled in the art will readily see that the precise 
steps of calculation and looping can be arranged in any 
desired or suitable way. An example of this is shown in 
Figs. 2a and 2b which show a second embodiment of the 
present invention. Steps in Figs. 2a and 2b which are the 

15 same as those already described in Figs, la and lb carry 
the same reference numerals and shall not be described 
again in detail. 

Instead of one loop S2, S3 and S4 for measuring Q^, Qrefj_ 
20 and determining x±, as in the case of Fig. 1, Fig. 2 has 
two loops S20, S3, S4 and S22, S23, S24. In other words, 
the measurement of and Qrefi and the determination of Xj_ 
is performed in two subsequent loops. On the other hand, 
the value R = w^ 2 is not determined in steps S9 and S10, as 
25 in Fig. 1, but only in one step S25, which precedes the 

decision in step S8. Therefore, as already mentioned, the 
present invention can basically be implemented in any 
desired way, as long as the sum 

30 R = W } 2 + X W j (6) 

7=2 

is determined in a recursive way by determining wj from the 
preceding value w-j_i. 

35 This is exemplified by a third embodiment of the present 
invention, a flow chart of which is shown in Fig. 6a and 
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5 6b. The third embodiment relates to a process that 

determines a deviation value TDEV. In Fig. 6a, step S100 
represents the determination of N values of xi, similarly 
to the above description of the first and second 
embodiments. Then, in step S101, n is assigned a value 

10 between 1 and the integer part of N/3 (referred to as 

Int(N/3)) . In the first two embodiments n was a given value 
with respect to the process for determining R, e.g. n was 
set by an external control means. As indicated in the 
present embodiment, the invention is by no means restricted 

15 thereto, and n can also be set in the process for 

determining R. It is clear that the setting of n is of no 
importance to the basic concept of the invention, namely 
the recursive determination of R, as n is only a parameter 
in said recursive determination process. 



20 



25 



In step S102 and the subsequent loop consisting of steps 
S103 to S105, similarly to steps S5 to. S8 in the previous 
embodiments, a first value of w is determined in accordance 
with 



n 



30 



Then, in step S106, a value of R is determined with the 
help of the above calculated value of w, namely R=w . If 
the value of n is equal to Int(N/3), then the calculation 
of R is complete, as step S107 will then branch to step 
S113 (see Fig. 6b) to calculate TDEV on the basis of the 
shown equation. On the other hand, if n is smaller than 
Int(N/3), then R will be calculated recursively in 
35 accordance with 
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W-3/7 + 1 9 

5 R = W 2+ W . Z (6) 

j = 2 

and 

wj = wj_i + acj_i + 3 n - 3xj_i +2n + 3xj_ 1+n - x-j-i (7) 

10 in steps S108 to S112. 

Finally, after R has been determined, TDEV is calculated in 
step SI 13 in accordance with the shown equation. 

15 A difference between the first two embodiments and the 
third embodiment is that in the third embodiment no 
vectorisation of w is necessary. 



20 



Figs. 3a and 3b show a flow-chart of a fourth embodiment of 
the present invention. Steps that are the same as in Fig. 1 
carry the same reference numerals and shall not be 
described again in detail. The first loop S2, S3 and S4, 
which is identical to the one shown in Fig. 1, determines 
error values Xi which characterize the variation between a 
25 quantity Q and a reference quantity Qref . The present 

embodiment is characterized by the definition of a help 
variable y k , as shown in the loop consisting of steps S51, 
S52 and S53 . As can be seen, in a first step S50, the loop 
variable k is set to 1, and then the help variable y k is 
30 defined in accordance with 



Yk = x k+2n " 2x n+k + x k 



(8) 



for values of k from 1 to N-2n. In the presently described 
35 embodiment, this is achieved by determining in step S52 if 
k is equal to N-2n, and if not, then incrementing k to k+1 
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in step SS3, after which the processing returns to step 
S51. These values of yjc are then stored. 

After all of the values y k for k = 1 to N-2n have been 
determined and stored, the processing enters a further loop 
referred to as S55, S56, S57 . In this loop/ which again 
runs through a loop variable referred to as k, the value w 2 
is determined in accordance with the equation (9) 



15 



20 



25 



30 



n 



(9) 

k=l 



35 



At first, in step S54, the loop variable k is set to 1 and 
the value w 1 is set to 0. Then vt 1 is replaced by v/i+Yk 
(Step S55) after which step S56 asks if k=n, and if not, 
then k is incremented to k + 1 in step S57, after which the 
processing returns to step S55. 

As in the case of in Fig. 2, the processing in Fig. 3 then 
continues with step S25, in which R is set equal to w^, 
after which step S81 determines if n is equal to the 
integer part of N/3. If not, then the loop formed by steps 
S60, S13, S14 and S15 is entered, which is comparable to 
the loops S12, S13, S14 and S15 known from the embodiments 
described in connection with Fig. 1 and 2, only that step 
S12 has been replaced by S60, in which the recursive 
determination of Wj is accomplished in accordance with 

wj = wj_i + Yj-l+n - Yj-1 (10) 

Although the above described embodiment contains the 
additional loop S51, S52 and S53 for determining the help 
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5 variable y^, the resulting calculation in the recursive 
equation of step S60 is simplified, which can be 
advantageous because it requires simpler hardware, 
especially if the method of the present invention is 
implemented in the form of circuitry containing adders, 

10 subtracters, etc. On the other hand, if the invention is 

implemented as software in an appropriate computing device, 
e.g. a CPU, then the embodiment described in connection 
with Fig. 1 is advantageous over that described in 
connection with Fig. 3, because the method embodied by the 

15 flow-chart of Fig. 1 contains one loop less than that of 
Fig. 3. 

As already mentioned in connection with Fig. 1, the value R 
determined by the methods described in connection with Fig. 
20 2 and 3 can also be further employed in any suitable and 
desired way, e.g. output to a control means that is 
controlled in accordance with the variation between Q and 
Qref, or processed further, e.g. for calculating a variance 

_ R 

or a deviation 



30 



DEV = 



R 



(12) 



The value R determined by the methods in accordance with 
the above described embodiments is equal to the double sum 
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J 



-i2 



L ' 



(3) 



described in the introduction. As already mentioned, the 
value n, which can be assigned any value between 1 and the 
integer part of N/3 describes the number of sampling 

10 intervals within one observation interval. This is 

explained in connection with Fig. 4, in which a number N of 
sampling intervals along the direction of a parameter t 
(for example time) are shown, each interval having a length 
of to- In tne specific example of Fig. 4, N = 12. The 

15 observation interval nx 0 relates to the triplets of values 
x associated with each sampling interval, i.e. Xi+2n' x i+n' 

*i- 



20 



25 



According to a preferred embodiment, the method of the 
present invention is employed such that all values of R 
relating to individual values of n, i.e. for all values of 
n from 1 to the integer part of N/3 are calculated. These 
values can then be processed further, e.g. for calculating 
the deviation DEV according to 



EEV = 



R 



6a 2 (Ar-3rc+l) 



(12) 



for each value of n from 1 to the integer part of N/3, in 
order to use these deviation values as appropriate 
30 characterization parameters. 



Such a measurement of all of the possible values R can 
simply be accomplished by introducing a further loop, as 
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5 e.g. done in Fig. 5, which is a modification of the 

embodiment shown in Fig. 1. All of the steps that are the 
same as in Fig. 1 carry the same reference numerals and 
shall not be described again. In addition to the steps of 
Fig. 1, the embodiment for calculating the values for n 

10 from 1 to the integer part of N/3 has the supplementary 
step S70 between steps S3 and S5, in which the loop 
parameter n is initially set to 1, and the supplementary 
steps S71, S72 and S73, where the value R calculated for a 
specific value n is stored in association with n in step 

15 S71, and steps S72 and S7 3 for the outer loop running over 
n, where S72 asks if n has reached the maximum value which 
is the integer part of N/3 and if not, then the processing 
proceeds to steps S73, in which n is incremented to n+1. 

20 After step S73, the loop returns to S5. As in the cases of 
the previous embodiments, the values R can be processed or 
output in any suitable or desired way. . 

It is clear that the embodiment described in connection 
25 with Figs. 5a and 5b can also be realized in other ways, 
for example also in connection with the embodiments 
described in Fig. 2 and 3. The modification of the method 
shown in Fig. 2 would consist in introducing the step S7 0 
shown in Fig. 5a between steps S23 and S5 and Fig. 2a, and 

30 adding steps S71, S72 and S73 at the end of the processing 
in Fig. 2b, as done in Fig. 5b. Similarly, the method 
described in connection with Fig. 3a would be modified by 
adding step S70 between steps S52 and S54 of Fig. 3a and 
adding steps S71, S72 and S73 at the end of the processing 

35 in Fig. 3b. This is readily understandable to a person 
skilled in the art. 

Fig. 7a and 7b show a modification of the embodiment of 
Fig. 6, where steps S201, S202 and S203 have been added in 
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20 



25 



30 



order to implement a calculation of TDEV (n x 0 ) for all 
values of n from 1 to the integer part of N/3. As can be 
seen, the loop for n is initialized in step S201 and loops 
back in step S202, and step S107 exits the loop when n is 
equal to the integer part of N/3. 



As already mentioned in connection with the embodiment of 
Fig. 1, the precise order and implementation of 
calculations and loopings in the embodiments of Figs. 2, 3, 
5, 6 and 7 is not important to the invention as long as the 
15 value R is calculated by recursively determining w in 
accordance with 

w . = wjfwj-!), i.e. from the previous value of w, and as a 
function of the difference values x, e.g. either by 



Wj = Wj_i + xj_i+3n " 3xj_i+2n + 3x j-l+n - (7) 

or by 

wj = wj_i + Yj-l+n ~ Yj-1 (10) • 

As can be seen in the embodiment of Fig. 5, the calculation 
of the values R for all values of n leads to one nested 
loop, in contrast to the conventional calculation which 
contains two nested loops, namely one for the inner sum 



£X (Xi+Zn-^+Xi) 



(3) 



and a second one for the outer loop running over n. Again, 
the number of calculations is reduced by a factor of N in 
35 the present invention, as the number or required 
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5 calculations in accordance with the embodiment shown in 
Fig. 5 is in the order of N 2 whereas the conventional 
calculation is in the order of N 3 . 

As a best mode of putting the invention to practice, the 
10 inventor presently considers the method described in 

connection with Fig. 7a and 7b, where the quantity Q is a 
clock signal T, Qref is a reference clock signal Tref, such 
that the values Xi are time error samples. The values R are 
calculated for all possible values of n, i.e. from n = 1 to 
15 n = the integer part of N/3, such that a complete set of 
values of the time deviation 



TDW{nT 0 ) = 



i 



l 



l 



6rt*(N-3rc+l) p 



>=j 



(2) 



20 can be calculated from the values R through 



EEV = 



i 



R 



(12) 



for all possible values of n. These values TDEV (nT 0 ) can 
25 then be used for characterization or supervision of 
reference timing signals. 



The precise implementation of the present invention in 
terms of hardware or software can be done in any suitable 
30 or desired way. For example, if a used software language 

has a provision such that if the running variable of a FOR- 
loop is outside of the range given by the upper and lower 
values of the FOR-loop, said FOR-loop is not entered, then 
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e.g. the decision step S81 and loop consisting of Sll, S14 
and S15 in Fig. 2b can be implemented by such a FOR-loop 
running from k=2 to N-3n+l. In other words, the decision 
step S81 would be implicitly contained in the above 
mentioned FOR-loop. 



Although the present invention has been described in terms 
of specific embodiments, it is by no means restricted 
thereto. Much rather, the scope of the invention is defined 
by the appended claims. Reference symbols in the claims 
15 serve the purpose of better understanding and do not 
restrict the scope. 
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5 [Claims] 

1. Method for measuring a value that characterizes a 
variation between a quantity (Q) and a reference 
quantity (Qref) associated with said quantity, both 
10 depending on a predetermined parameter (t), 

comprising: 

- measuring N samples (Qi, Qrefi) of said quantity (Q) 
and said reference quantity (Qref) at equally spaced 
15 intervals (To) of said parameter (t), N being an 

integer with N > 3, 



- calculating and storing N difference values xi 
between each sample of said quantity (Qi) and each 

20 sample of said reference quantity (Qrefi) at each of 

said N sample points, the index i being an integer 
running from 1 to N, 

- calculating a value wi as a function of said 
25 difference values xi, 

- calculating a value R representative of the 
variation between said quantity and said reference 
quantity as 

30 

N-3?i + \ 2 
R = w 2 + £ Wj , 

7 = 2 

where n is an integer from the set of 1 to the integer 
part of N/3, and each value Wj with j>l is calculated 
35 recursively from the previous value Wj_ x and as a 
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5 function of said difference values x it if n is smaller 

than the integer part of N/3, and 

- calculating said value R as R = wi 2 if n is equal to 
the integer part of N/3. 

10 

2. Method according to claim 1, characterized by 

- calculating said value w x according to the equation 



15 



k = \ 



and 

20 - calculating said each value w-j with j>l from the 

previous value Wj_i in accordance with 

25 3. Method according to claim 1, characterized by 
- calculating and storing values y k defined as 
v, —x, „ —2x , , + X, , 

Sk k+2n n+k k 



30 



for k= 1 to N-2n, 

- calculating said value w x according to the equation 
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n 

w, = 2>* > and 



- calculating said each value Wj with j>l from the 
previous value Wj-j in accordance with 



10 W^Yl+V""^' 



The method of one of claims 1 to 3, characterized in 
that a variance value VAR is calculated as 



R 

VAR 



6n 2 (N-3n+\) 

5. The method of claim 4, characterized in that said 
variance value is calculated for each value of n from 
1 to the integer part of N/3. 

6. The method of one of claims 1 to 3, characterized in 
that a deviation value DEV is calculated as 



DEV ■ R 



6n 2 (N-3n+l) 

7. The method of claim 6, characterized in that said 

deviation value is calculated for each value of n from 
1 to the integer part of N/3. 

30 8. The method of one of claims 1 to 7, characterized in 
that said parameter on which said quantity and said 
reference quantity depend is time. 
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9. The method of claim 8, characterized in that said 
quantity is a timing signal and said reference 
quantity is a reference timing signal. 

10. The method of one of claims 1 to 7, characterized in 
that said parameter on which said quantity and said 
reference quantity depend is location. 

11. The method of claim 10, characterized in that said 
quantity is a surface profile level and said reference 
quantity is a reference profile level. 

12. A device for measuring a value that characterizes a 
variation between a quantity (Q) and a reference 
quantity (Qref) associated with said quantity, both 
depending on a predetermined parameter (t) , 
comprising: 

- a measuring means for measuring N samples (Qi, 
Qrefi) of said quantity (Q) and said reference 
quantity (Qref) at equally spaced intervals (to) of 
said parameter (t) , N being an integer with N ;> 3, 

- a calculating and storing means for calculating and 
storing N difference values Xi between each sample of 
said quantity (Qi) and each sample of said reference 
quantity (Qrefi) at each of said N sample points, the 
index i being an integer running from 1 to N, 

- a calculating means for calculating a value w x as a 
function of said difference values x if for calculating 
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a value R representative of the variation between said 
quantity and said reference quantity as 



where n is an integer from the set of 1 to the integer 
part of N/3, and each value with j>l is calculated 
recursively from the previous value Wj_i and as a 
function of said difference values xi, if n is smaller 
than the integer part of N/3, and for calculating said 
value R as R = wi 2 if n is equal to the integer part 
of N/3. 

13. Device according to claim 12, characterized by said 

calculating means calculating said value w x according 
to the equation 



R = w 2 + £ Wj , 

j = 2 



n 




k = l 



and calculating said each value Wj with j>l from 
previous value Wj_i in accordance with 



W . = W 
J 



+x 



7-l+3» 



-3x 



7-1+2/j 



+ 3x 



j-l+n 



— X 



7-1 * 



Device according to claim 12, characterized by said 
calculating means calculating and storing values yk 
defined as 
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v, = x, „ —2x . +X. , 

*k k+2n n+k k 

for k= 1 to N-2n, calculating said value wi according 
to the equation 

n 

w, = Hyk ' and 

k = 1 

calculating said each value Wj with j>l from the 
previous value Wj_i in accordance with 
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Fig. 3a 
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Fig. 4 
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